PARAMETER MISMATCHES, VARIABLE DELAY TIMES AND SYNCHRONIZATION IN 

TIME-DELAYED SYSTEMS 
E. M. Shahverdiev 1 , R.A.Nuriev 
Institute of Physics, 370143 Baku, Azerbaijan 
R.H.Hashimov 
Azerbaijan Technical University, 370073 Baku, Azerbaijan 

K. A. Shore 

School of Informatics, University of Wales, Bangor, Dean Street, Bangor, LL57 1UT, Wales, UK 

ABSTRACT 

We investigate synchronization between two unidirectionally linearly coupled chaotic non-identical 
time-delayed systems and show that parameter mismatches are of crucial importance to achieve 
synchronization. We establish that independent of the relation between the delay time in the 
coupled systems and the coupling delay time, only retarded synchronization with the coupling 
delay time is obtained. We show that with parameter mismatch or without it neither complete 
nor anticipating synchronization occurs. We derive existence and stability conditions for the re- 
tarded synchronization manifold. We demonstrate our approach using examples of the Ikeda and 
Mackey Glass models. Also for the first time we investigate chaos synchronization in time-delayed 
systems with variable delay time and find both existence and sufficient stability conditions for the 
retarded synchronization manifold with the coupling- delay lag time. Also for the first time we 
consider synchronization between two unidirectionally coupled chaotic multi-feedback Ikeda sys- 
tems and derive existence and stability conditions for the different anticipating, lag, and complete 
synchronization regimes. 
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1. INTRODUCTION 

Seminal papers on chaos synchronization [1] have stimulated a wide range of research activity 
especially extensively in lasers, electronic circuits, chemical and biological systems [2]. Possible 
application areas of chaos synchronization are in secure communications, optimization of nonlinear 
system performance, modeling brain activity and pattern recognition phenomena [2]. 

There are different types of sychronization in interacting chaotic systems. Complete, gen- 
eralized, phase, lag and anticipating synchronizations of chaotic oscillators have been described 
theoretically and observed experimentally. Complete synchronization implies coincidence of states 



Regular Associate with the Abdus Salam ICTP corresponding author:shahverdiev@physics.ab.az 

1 



of interacting systems, y(t) = x(t) [1]. Generalized synchronization is defined as the presence of 
some functional relation between the states of response and drive, i.e. y(t) = F(x(t)) [3]. Phase 
synchronization means entrainment of phases of chaotic oscillators, n$ x — m& y = const, (n and m 
are integers) whereas their amplitudes remain chaotic and uncorrelated [4]. Lag synchronization 
for the first time was introduced by Rosenblum et al. [5] under certain approximations in study- 
ing synchronization between bi-directionally coupled systems described by the ordinary differential 
equations (no intrinsic delay terms) with parameter mismatches:y(t) ~ x T (t) = x(t — r) with pos- 
itive t. Anticipating synchronization [6-8] also appears as a coincidence of shifted-in-time states 
of two coupled systems, but in this case the driven system anticipates the driver, y(t) = x(t + r) 
or x = y T ,r > 0. An experimental observation of anticipating synchronization in external cavity 
laser diodes [9] has been reported recently, see also [10] for the theoretical interpretation of the 
experimental results. The concept of inverse anticipating synchronization x = —y T is introduced 
in [11]. 

Due to finite signal transmission times, switching speeds and memory effects time-delayed 
systems are ubiquitous in nature, technology and society [12]. Therefore the study of synchro- 
nization phenomena in such systems is of high practical importance. Time-delayed systems are 
also interesting because the dimension of their chaotic dynamics can be made arbitrarily large by 
increasing their delay time. From this point of view these systems are especially appealing for 
secure communication schemes [13]. 

Role of parameter mismatches in synchronization phenomena is quite versatile. In certain cases 
parameter mismatches are detrimental to the synchronization quality: in the case of small param- 
eter mismatches the synchronization error does not decay to zero with time, but can show small 
fluctuations about zero or even a non-zero mean value; larger values of parameter mismatches can 
result in the loss of synchronization [8,14]. In some cases parameter mismatches change the time 
shift between the synchronized systems [15]. In certain cases their presence is necessary for syn- 
chronization. We reiterate that the crucial role of parameter mismatches for lag synchronization 
between bi-directionally coupled systems was first studied in [5] by Rosenblum et al. As such, lag 
synchronization cannot be observed if two oscillators are completely identical, see e.g. [16] and 
references therein. 

Multi-feedback and multi-delay systems are ubiquitous in nature and technology. Prominent 
examples can be found in biological and biomedical systems, laser physics, integrated communi- 
cations [12]. In laser physics such a situation arises in lasers subject to two or more optical or 
elctro-optical feedback. Second optical feedback could be useful to stabilize laser intensity [17]. 
Chaotic behaviour of laser systems with two optical feedback mechanism is studied in recent works 
[18]. To the best of our knowledge chaos synchronization between the multi-feedback systems is to 
be investigated yet. Having in mind enormous application implications of chaos synchronization 
e.g. in secure communication, investigation of synchronization regimes (lag, complete, anticipat- 
ing etc.) in multi-feedback systems is of immense importance. 

In this paper we investigate synchronization between the two unidirectionally coupled chaotic 
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non-identical time-delayed systems having a fairly general form of coupling and show for the first 
time that parameter mismatches are, in fact, of crucial importance for achieving synchroniza- 
tion. We show that independent of the relation between the delay time in the coupled systems 
and the coupling delay time, only retarded (lag) synchronization is obtained. ( Usually for lag 
synchronization between the unidirecitionally coupled time-delayed systems the term retarded 
synchronization is prefered [8].) In this case the lag time is the coupling delay time. We consider 
both constant and variable feedback delay times. We demonstrate our approach using examples 
of the Ikeda and Mackey Glass models. 

In the paper also for the first time we investigate synchronization between two unidirectionally 
coupled chaotic multi-feedback Ikeda systems and find both existence and stability conditions for 
different synchronization regimes. 



2. GENERAL THEORY 
Consider a situation where a time-delayed chaotic master (driver) system 



clx 

— = -a 1 x + k 1 f(x T1 ), (1) 

drives a non-identical slave (response) system 

dy 

— = -a 2 y + k 2 f(y T1 ) + k 3 x T2 , (2) 

where x and y are dynamical variables; f(x) is different iable nonlinear function; a.\ and a 2 are 
relaxation coefficients for the driving and driven dynamical variables, respectively:throughout the 
paper we assume that ot,\ = a — 5 and a 2 = a + 5, 5 determines the mismatch of relaxation 
coefficents; T\ is the feedback delay time in the coupled systems;r2 is the coupling delay time 
between the systems. k\ and k 2 are the feedback rates for the master and the response systems, 
respectively;/^ is the linear coupling rate between the driver and the response system. 
Now we will show that chaotic systems (1) and (2) can be synchronized on the retarded synchro- 
nization manifold with the lag time r 2 . 

y = x T2 . (3) 

We denote the error signal by A = x T2 — y. Then from systems (1) and (2) we find the following 
error dynamics: ^ = — a 2 A + {25 — k 3 )x T2 + kif(x Tl+T2 ) — k 2 f(y Tl ). Thus under conditions 

25 = k 3 , h = k 2 , (4) 

the error dynanics can be written as: 

^ = -a 2 A + A; 1 A T1 / / (x ri+T2 ). (5) 
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It is obvious that A = is a solution of system (5). To study the sufficient stability condition 
for the retarded synchronization manifold y = x T2 one can use a Krasovskii-Lyapunov functional 
approach [12, 19]. 

The sufficient stability condition for the trivial solution A = of eq.(5) can be found by investi- 
gating the positively defined Krasovskii-Lyapunov fuctional 

V(t) = iA 2 + /i J^A^t + tJdh, (6) 

where \x > is an arbitrary positive parameter. According to [12,19], the solution A = is stable, 
if the derivative of the functional (6) along the trajectory of equation ^ = — r(t)A — s(t)A T is 
negative. In general this negativity condition is of the form: 4(r — fij/j, > s 2 and r > fi > 0. As 
the value of fj, that will allow s 2 as large as possible is fj, — |, the asymptotic stability condition 
for A = can be written as 

r 2 > s 2 , (7) 

which is equivalent to r > \s\. This result is valid for both constant and time-dependent coefficients 
r and s (in the latter case r(t) and s(t) should be bounded continuous functions [12]). Thus we 
obtain that 

a 2 > \k-tf {x Tl+T2 )\ (8) 

is the sufficient stability condition for retarded synchronization manifold (3). The condition (4) is 
the existence condition of retarded synchronization between the unidirectionally coupled systems 
(1) and (2). 

Thus we find that under certain conditions systems (1) and (2) admit the retarded chaos synchro- 
nization manifold y = x T2 only under parameter mismatch ie ol\ ^ ol 2 - We also notice that without 
the parameter mismatch, i.e. a.\ = a 2 = a neither y = x T2 - T1 nor y = x Tl - T2 is the synchronization 
manifold. We also emphasize that, in general for both cci = a 2 and a.\ 7^ a 2 systems (1) and (2) 
admits neither complete nor anticipating chaos synchronization. 

So far we have considered the case of constant feedback delay time t±. It is of immense interest 
to study chaos synchronization in time-delayed systems with variable feedback delay time. Basic 
interest is driven by the fact that so far there are no reported research on this particular subject in 
the literature. Practical interest is motivated by the appreciation that time-delayed systems with 
variable delay times are more realistic. As an example one can refer to the biological biorhythms, 
where the capacity of assimilation of nutrients by an organism varies cyclicly during the day [20] . 
Now we will try to find both the existence and stability conditions for the synchronization man- 
ifold (3) in the case of variable feedback delay times. It is straightforward to establish that the 
analog of the error dynamics equation in the case of variable delay time T\{t) is of the form: 

^ = -a 2 A + A; 1 A ri(t) / / (x Tl(t)+r2 ). (9) 
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Again, as in the case of constant feedback delay times equation (9) is obtained from the studying 
the coevolution of eqs.(l) and (2) along the manifold (3). Analysis of the error dynamics shows 
that the existence conditions (4) hold for the variable delay cases. Next let us find the sufficient 
stability condition for system (9). According to [12] for that purpose one can still use the functional 
(6). Namely as presented in [12], when r = r(t) is continuously differentiable and bounded, the 
solution A = to ^ = — r(t) A — s(t) A T ( t ) is uniformly asymptotically stable, if a(t) > pL > and 
(2r(t) — yu) (1 — %)a* > s2 (^) uniformly in t. Applying the same procedure as in the case of constant 
feedback delay time, we can find the value of fi that will allow s 2 to be as large as possible: fi — r. 
Thus we find that the sufficent stability condition for the A = solution of time delay equation 
with time dependent coefficients ^ = —r(t)A — s(t)A r (t) is: 

r\t)(l-^)>s 2 (t). (10) 

Notice that for the constant delay time cases the inequality (10) is reduced to the well-known 
sufficent stability condition r > \s\. 

As in our case r = ct 2 and s = —kif'(x T1 ^ +T2 ) then the sufficent stability condition for synchro- 
nization manifold (3) for the time-delayed equations (1) and (2) with time dependent feedback 
delay T\ can be written as: 

^ (1 "lf )>(fcl/ ' ( ^ (t)+T2))2 (U) 
3.1. EXAMPLE 1:THE IKEDA MODEL 

In this subsection we demonstrate our general theory using the example of the Ikeda model. This 
investigation is of considerable practical importance, as the equations of the class B lasers with 
feedback (typical representatives of class B are solid-state, semiconductor, and low pressure C0 2 
lasers [21]) can be reduced to an equation of the Ikeda type [22]. Consider synchronization between 
the Ikeda systems [6], 

dx 
~dt 
dy 

— = -a 2 y - psiny Tl + Kx T2 . (12) 

The Ikeda model was introduced to describe the dynamics of an optical bistable resonator and is 
well-known for delay-induced chaotic behavior [23]. Physically x is the phase lag of the electric 
field across the resonator; a is the relaxation coefficient; (3 is the laser intensity injected into the 
system. T\ is the round trip time of the light in the resonator or feedback delay time in the coupled 
systems;r 2 is the coupling delay time between systems x and y. 

First we consider the case of constant feedback delay time and show that y = x T2 is the retarded 



-a.\x — /3sin x T1 , 
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synchronization manifold, if the parameter mismatch a 2 — &i = 2<5 is equal to the coupling rate 
K. This can be seen by the dynamics of the error A = x T2 — y: 



dA 
~dt 



(a + 5) A + (25 — K)x T2 — (3 cos x. 



A 



ri ■ 



(13) 



(As in this example under study we choose feedback rates ((3) equal for both the driver and driven 
systems, the second of the existence conditions in (4) becomes redundant.) The sufficient stability 
condition for the retarded synchronization manifold y = x T2 can be written as:a + 8 — a 2 > |/3|. 
Thus, as in case of general approach, we find that the retarded chaos synchronization manifold 
y = x T2 occurs only under parameter mismatch ie ct\ ^ a 2 . By analyzing the corresponding 
error dynamics one can also establish that without the parameter mismatch, i.e. ot\ = a 2 = a 
neither y = x T2 - T1 nor y = x T1 _ T2 is the synchronization manifold. We also emphasize that for 
both cci = a 2 and a± ^ a 2 system (12) admits neither complete (we notice that for special 
case of r 2 = y = x T2 is the complete synchronization manifold, which exists if ol\ ^ a 2 ) nor 
anticipating chaos synchronization. We emphasize that this result is due to the linear coupling 
between the synchronized systems. The importance of the role of the form of coupling between 
the synchronized systems is underlined in [6,24]. In the case of nonlinear (sinusoidal) coupling for 
identical drive and response Ikeda systems, depending on the relation between the feedback delay 
time and the coupling delay time retarded, complete or anticipating synchronization can occur, 
see, e.g. [25] and references therein. 

According to estimations the parameters values ct\ = 5, a 2 = 25, (3 = 20, K = 20 and T\ = 1, 
t 2 = 2 (or Ti = 3, t 2 = 1) satisfy both existence and stability conditions for the retarded chaos 
synchronizartion manifold y = x T2 for the coupled systems (12). 

Next we consider the case of time dependent delay time Ti(t). First we notice that as in 
the case of time-independent delay times 26 = K is the condition of existence for the y = x T2 
synchronization manifold. Next applying the general formula (11) derived earlier in the paper we 
write the sufficent stability condition for the synchronization manifold y = x T2 in the following 
form: 



where To is the zero frequency component;r a is the amplitude; ^ is the frequency of the modulation. 
Then for the concrete form of variable delay time (15) the sufficient stability condition (14) can 
be writen as: 




(14) 



As an example consider the following sinusoidal form of the variable delay time : 



Ti(t) =r + T a s'm(ujt), 



(15) 



a|(l — r a uj cos(c<jt)) > /3 2 . 
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(16) 



3. 2. EXAMPLE 2:THE MACKEY GLASS MODEL 



In this subsection we demonstrate our approach using the example of the Mackey Glass model. 
The Mackey Glass model has been introduced as a model of blood generation for patients with 
leukemia and nowadays is very popular in chaos theory [19]. 
Consider synchronization between the Mackey Glass systems 

— — —OL\X + K\- — ; — , 



dt 1 + x 



Tl 



dy a 2 x T 
— = -a 2 y + k 2 —— h 
dt 1 + x° 



-a 2 y + k 2 - r —^- + hx T2 . (17) 

The dynamical variable in the Mackey Glass model is the concentration of the mature cells in 
blood at time t and the delay time is the time between the initiation of cellular production in the 
bone marrow and the release of mature cells into the blood [20]. 

Again by investigating the corresponding error dynamics we can show that y = x T2 is the retarded 
synchronization manifold, if the parameter mismatch a 2 — a.\ = 25 is equal to the coupling 
rate k% and k\a\ = k 2 a 2 . We notice that here we can allow for parameter mismatches for a, 
and thus have more flexibility to achieve synchronization. With these existence conditions, the 
sufficient stability condition for the retarded synchronization manifold y = x T2 can be written as: 
a 2 > \kiaif'(x Tl+T2 )\, with f(x T ) = 

For analytical estimation of a 2 we take into account that the absolute maximum of the function 
\f'(x T )\ is obtained at x T = (frf)^ and is equal to ^ 6 ~^ [19]. Thus we arrive at the following 
sufficient stability condition for the synchronization manifold y = x T2 for the coupled systems (17) 

a 2 > Mi Ab ' • (18) 

Again we would like to underline that only retarded synchronization occurs notwithstanding the 
relation between the feedback delay time and coupling delay time; moreover for both a.\ = a 2 and 

7^ ct 2 coupled systems (17) admit neither complete nor anticipating chaos synchronization. 
According to estimations the parameters values cci = 0.1, a 2 — 5, k±ai = 2,6 = 10, k 3 = 4.9,Ti = 10 
and t 2 = 20 (or T\ = 50 and r 2 = 20, T\ = r 2 = 50) satisfy both existence and stability conditions 
for the retarded chaos synchronizartion manifold y = x T2 for the coupled systems (17). 

4. SYNCHRONIZATION BETWEEN THE MULTI-FEEDBACK IKEDA SYSTEMS 
Consider synchronization between the multi-feedback Ikeda systems, 

dx 

— = —ax + mi s m x n + "^2 sin x T2 , (19) 
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dy 

— = —ay + m 3 sin y T1 + m 4 sin y T2 + K sin x TZ , (20) 

CLL 

with positive ai, 2 and —7711,2,3,4. 

As mentioned above physically x is the phase lag of the electric field across the resonator; a is 
the relaxation coefficient for the driving x and driven y dynamical variables;— 777,1,2 and — 7773,4 are 
the laser intensities injected into the driving and driven systems,respectively. 71,2 are the feedback 
delay times in the coupled systems;r3 is the coupling delay time between systems x and y;K is 
the coupling rate between the driver x and the response system y. 
First we will show that systems (19) and (20) can be synchronized on the manifold: 

y = x 7i - T1 . (21) 

(For r 3 > Ti,t 3 = Ti, and t 3 < T\ (21) is the retarded, complete, and anticipating synchronization 
manifold, respectively.) We denote the error signal by A = x T:i ^ Tl — y. Then from systems (19) 
and (20)we find the following error dynamics A = x T3 _ Tl — y'^rr = —cxA + ((mi — K) sinx T3 — 
m 3 sin y Tl ) + ?77 2 sin x T2+T3 - Tl — 7774 sin y T2 Thus under conditions 

777i — K = 7773, 777 2 = 7774 (22) 

the error dynamics can be written as: 

— — = -aA + 777 3 A Tl cosx T3 + 777 2 A T2 cosx T2+T3 _ ri . (23) 

(JjL 

It is obvious that A = is the solution of system (23). To study the stability of the synchronization 
manifold y = x T3 _ T1 one can again use a Krasovskii-Lyapunov functional approach. According to 
[12], the sufficient stability condition for the trivial solution A = of time-delayed equation 
^ = -r(t)A + Sl (t)A Tl + s 2 (t)A T2 is: r(t) > \ Sl (t)\ + \s 2 (t)\. 

Thus we obtain that the sufficient stability condition for the synchronization manifold y = x T3 _ Tl 
(21) can be written as: 

a > 1 777 3 1 + 1 7772 1 . (24) 

Conditions (22) are the existence conditions for the synchronization manifold (21) between the 
unidirectionally coupled multi-feedback systems (19) and (20). 

Analogously we find that y = x T3 _ T2 is the synchronization manifold between systems (19) and (20) 
with corresponding existence 777 2 — = 7774 and 777,1 = 7773 and stability conditions a > \m 3 \ + \ m^\. 

One can easily generalize the previous rezults to 77-tuple feedback Ikeda systems. Indeed con- 
sider synchronization between the following Ikeda models: 

dx . . . 

— = —ax + 777i x sin x T1 + 777 2x sin x T2 + ■ ■ ■ + m nx sin x Tn , (25) 



— = -ay + m ly sin y T1 + m 2y sin y T2 + m ny sin y Tn + k sin x Tk , (26) 

Then the existence and sufficent stability conditions e.g. for the synchronization manifold y = 
x Tk _ T1 are:m lx — k = mi y ,m nx = m ny and a > \rrii y \ + \m 2y \ + • • • + |m n?/ |, respectively. For 
synchronization manifold y = x Tk _ T2 , m 2x — k = m 2y and m nx = m ny are the existence conditions, 
and a > \mi y \ + \m 2y \ + • • • + \m ny \ is the sufficient stability condition. 



5. CONCLUSIONS 

In this paper we have studied the relation between parameter mismatches and synchronization 
in a certain class of unidirectionally linearly coupled time-delayed systems and have shown for 
the first time that parameter mismatches are of crucial importance for achieving synchronization. 
We have showed that independent of the relation between the feedback delay time in the coupled 
systems and the coupling delay time, only retarded (lag) synchronization with coupling delay lag 
time is obtained. We have established that either with parameter mismatch or without it neither 
complete nor anticipating chaos synchronization occurs. We have demonstrated our approach 
using the Ikeda and Mackey Glass models. We mention that, for example in the case of nonlinear 
(sinusoidal) coupling for identical drive and response Ikeda systems, depending on the relation 
between the feedback delay time and the coupling delay time retarded, complete or anticipating 
synchronization can occur [25]. These results are of significant interest in the context of relation- 
ship between parameter mismatches, coupling forms and synchronization. Indeed, having in mind 
possible practical applications of anticipating chaos synchronization [6] in secure communications 
(anticipation of the future states of the transmitter (master laser) at the receiver (slave laser) 
allows more time to decode the message ), in the control of delay-induced instabilites in a wide 
range of non-linear systems, for the understanding of natural information processing choosing the 
"appropriate" parameters' mismatches and coupling forms certain types of synchronization can 
be switched off/on. We have also for the first time investigated chaos synchronization in variable 
delay time systems and found both existence and sufficient stability conditions for the retarded 
synchronization manifold with the coupling-delay lag time. 

Also for the first time we have investigated synchronization between two unidirectionally coupled 
multi-feedback systems. These findings are of considerable interest in the context of synchro- 
nization between the stabilized laser systems (arrays) which hold great promise for space com- 
munication applications, where compact sources with high intensities are required. Additionally, 
synchronization between the multi-feedback systems can provide more flexibility in practical ap- 
plications. 
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